function jump=q(x)
jump=zeros(size(x));
jump(1,1)=2*(3*x(1,1)^2+2*x(1,1)*x(1,2)+x(1,2)^2)/sqrt(x(1,1)^2+x(1,2)^2);
jump(1,2)=2*(3*x(1,2)^2+2*x(1,1)*x(1,2)+x(1,1)^2)/sqrt(x(1,1)^2+x(1,2)^2);
